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ABSTRACT 

We present the updated version of the code used to compute stellar evolutionary tracks in 
Padova. It is the result of a thorough revision of the major input physics, together with the 
inclusion of the pre-main sequence phase, not present in our previous releases of stellar mod- 
els. Another innovative aspect is the possibility of promptly generating accurate opacity tables 
fully consistent with any selected initial chemical composition, by coupling the OPAL opac- 
ity data at high temperat ures to the molecular opacities computed with our ifiSOPUS code 
dMarigo & Aringerl2009l) . In this work we present extended sets of stellar evolutionary models 
for various initial chemical compositions, while other sets with different metallicities and/or 
different distributions of heavy elements are being computed. For the present rel ease of mod- 
els we adopt the solar distribution of heavy elements from the recent revision by Caffau et alj 
d201 lb . corresponding to a Sun's metallicity Z 0.0152. From all computed sets of stellar 
tracks, we also derive isochrones in several photometric systems. The aim is to provide the 
community with the basic tools to model star clusters and galaxies by means of population 
synthesis techniques. 
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1 INTRODUCTION 

In this paper we briefly describe our new stellar evolution code 
PARSEC: the PAdova & TRieste Stellar Evolution Code. This is 
the result of a thorough revision and update of the stellar evolution 
code used in Padova to compute sets of stellar evolu tionary tracks 
that a r e widely used by the astronomical community i Bressan et al.1 
19931: iBertellietal . 1 1 1994 iGirardi et alj|2000l : iMarigo et alj|200ll : 
Bertelli et al . 2008, 20q1). 



As we describe next, the most important changes include the 
updating of the major input physics (equation of state, opacities, 
nuclear reaction rates), the implementation of microscopic diffu- 
sion, and the extension of the evolutionary calculations to the pre- 
main sequence phase (PMS). 

We have devoted much care to ensure consistency between the 
adopted chemical mixture and the physical ingredients in PARSEC, 
i.e. opacities and equation of state, as detailed in Sects. l2~2l and !2.3l 
respectively. In this respect a distinctive prerogative of our stella r 
models is the use of the jESOPUS tool dMarigo & Aringerll2009r) . 
that enables us to compute tables of accurate low-J opacities for 
any specified mixture of chemical elements from H to U, just start- 



ing from the monochromatic absorption coefficients of all opacity 
sources under consideration (e.g. line lists of several molecules). 
We have put the most possible effort to keep the same abundance 
flexibility also when dealing with the equation of state, the high- 
temperature radiative opacities and the conductive opacities, for 
which we have employed other public codes. This represents a key 
aspect since we are now able to predict readily the evolution of stars 
for any chemical pattern of interest (varying CNO abundances, dif- 
ferent degrees of enhancement/reduction in a-elements, C-N, Ne-0 
and Mg-Al abundance anti-correlations suitable for Galactic glob- 
ular clusters, etc.). 

In this paper we present new sets of evolutionary tracks for 
the initial chemical compositions in the range from Z = 0.0005 
to Z = 0.07. The helium content Y is assumed to increase with Z. 
Other metallicities and arbitrary distributions of heavy elements are 
being considered. The range of initial masses ranges from 0. 1 M Q 
to 12 M , and the evolutionary phases extend from the PMS till 
either the onset thermally-pulsing AGB phase (TP-AGB) or carbon 
ignition. From all these tracks, we also derive isochrones in several 
photometric systems. 

The plan of this paper is as follows. Section|2]presents the in- 
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Table 1. Solar abundances of a few elements adopt ed in this work, follow- 
ing the values recommended by Cattau et alj 1201 ll) and references therein. 
Abundances are expressed with the standard notation A (Y) = \og(ny /«h) + 
12. Fo r all other species we adopt the compilation of IGrevesse & Sauvaj 
<1998t) . 



Element 


Abundance A(Y) 


Reference 


Li 


1.03 


Caffauetal. (2011) 


C 


8.50 


Caffauetal. (2010) 


N 


7.86 


Caffau et al. (2009) 





8.76 


Caffau et al. (2008) 


P 


5.46 


Caffau et al. (2007) 


s 


7.16 


Caffau & Ludwia (2007) 


K 


5.11 


Caffauetal. (2011) 


Fe 


7.52 


Caffau et al. (2011) 


Eu 


0.52 


Mucciarelli et al. (2008) 


Hf 


0.87 


Caffau et al. (2008) 


Os 


1.36 


Caffauetal. (2011) 


Th 


0.08 


Caffau et al. (2008) 



put physics of the models which includes the new solar distribution 
of elements; Section[3]deals with the calibration of important model 
parameters with the Solar model; Section|4]briefly describes other 
adopted element distributions and global metallicities; Section [5] 
introduces the new stellar tracks and discuss their main character- 
istics and Section[6]describes the corresponding isochrones. 



2 INPUT PHYSICS 

2.1 The solar distribution of heavy elements 

Before discussing any other relevant input physics we specify the 
solar distribution of heavy elements adopted in this paper. For each 
element heavier than 4 He, we must assign its fractional abundance 
relative to the total solar metallicity, i.e. X, /Z . This will be the 
reference distribution, with respect to which other non-scaled solar 
mixtures will be considered. 

The reference solar distribution of metals consists of 90 chem- 
ical elem ent j^from Li to U, with abu ndances taken from the compi- 
lation by IGrevesse & Sauvall dl998l) . except for a subset of species 
for which we adopt the recomme nded values according to the lat- 
est revision bv lCaffauet ai]j2fJT 11) and references therein. The solar 
abundances of the recently revised elements are listed in Table [fl 

According to this abundance compilation, the present-day 
Sun's metallicity is Z Q = 0.01524, tha t can be compared to other 
recent estimate s, e.g. Z n = 0.0141 of lLodders et alj d2009l) . and 
Z = 0.0134 of lAsplund et~ai] d2009l) . 

In addition to the solar one, other chemical distributions have 
been considered in this work, which should be representative of 
galaxies with specific chemical evolution histories such as, for in- 
stance, the Magellanic Clouds, massive early-type galaxies or other 
systems with ^-enhanced or fr-depleted patterns. 



2.2 Opacities 

Following a standard procedure, in PARSEC we describe the ab- 
sorption properties of matter in the gas phase by means of pre- 
computed, static tables of Rosseland mean opacities, K R (p, T), 
which are suitably arranged to encompass a region of the 
temperature-density (T - p) plane wide enough to cover all values 
met across the stellar structure during the evolutior0. 

Beside the state variables, T and p, opacities depend on the 
chemical composition of the gas, which is commonly specified by 
a set of abundance parameters, i.e.: the total metallicity Z, the hy- 
drogen abundance X, and the distributions [Xj/Z] of the heavy ele- 
ments in the mixture. The latter depend on the specific case under 
consideration, e.g. scaled-solar mixtures with [Xj/Z] = {X,-/Z} , or 
others [Xj/Z] derived from various constraints such as the enhance- 
ment/depletion of a-elements (expressed by the ratio [a/Fe]), or 
the over-abundances in primary C and O necessary to describe the 
He-burning regions. For the computation of the opacity tables we 
employ different programs as detailed below. 

In the high-temperature regime, 4.2 < log(T/K) < 8.7, 
we adopt the opacity tab les provided by the Opa city Project At 
Livermore (OPAL) team dlglesias & Rogersll 199(1 and references 
therein). We use the OPAL interactive web maslQ specifying the 
number fractions of 19 heavy elements (C, N, O, Ne, Na, Mg, Al, 
Si, P, S, CI, Ar, K, Ca, Ti, Cr, Mn, Fe, and Ni), as implied by our 
set of Xj/Z. 

In the low-temperature regime, 3.2 < log(r/K ) ^ 4.1, we 
employ the .ESOPUS^I tool dMarigo & Aringerll2009|) to generate 
opacity tables for any specified set of chemical abundances for 
92 elements from H to U. ^ESOPUS solves the equation of state 
of matter in the gas phase for ss 800 chemical species consisting 
of almost 300 atoms (neutral atoms and ions up to 5 th ionization 
stage) and 500 molecular species. The gas opacities account for 
many continuum and discrete sources, including atomic opacities, 
molecular absorption, and collision-induced absorption. 

In the transition interval 4.0 < log(T/K) < 4. 1, a linear inter- 
polation between the opacities derived from the OPAL and jESO- 
PUS is adopted. We remind the reader that both opacities sources 
provide v alues in good agreement in this temperature interval. As 
shown bv lMarigo & Aringej d2009l see their figure 7) the logarith- 
mic differences between OPAL and ^ESOPUS opacities are lower 
than 0.05 dex for most cases. 

Conductive opacities are included following lltoh et al.|j2008h . 
In the computation, for any specified chemical mixture, the total 
thermal conductivity accounts for the contribution of 11 atomic 
species (»H, 4 He, 12 C, 14 N, 16 0, 20 Ne, 24 Mg, 28 Si, 32 S, 40 Ca, and 
56 Fe), each weighted by the corresponding abundance (by number). 
To this aim we have implemented in our code the fortran routine 
kindly made available by Itoh (private communication). 

In practice, given the total reference metallicity Z and the dis- 
tribution of heavy elements {Xj/Z}, we construct two sets of opacity 
tables, to which we simply refer to as "H-rich" and "H-free" opaci- 
ties. The former set comprises N x opacity tables, where N x denotes 
the number (typically 10) of hydrogen abundance values, ranging 
from X = to X = 1-Z. The latter set is characterised by X = 0, 



A few elements (Po, At, Rn, Fr, Ra, Ac, and Pa) are assigned negligible 
abundances. 

2 Some of the elements listed in Table[T](e.g. Eu, Hf, Th) may be irrelevant 
for the evolutionary calculations presented in this paper, however they will 
become relevant in the context of the TP-AGB tracks, to be described in 
subsequent papers. 



3 Under the assumption of local thermodynamical equilibrium and in the 
limit of diffusion approximation, the solution of the radiation transfer equa- 
tion simplifies in such a way that the flux of radiation F as a function of 
radius r is expressed as F(r) oc 1/kr. 

4 http://opalopacity.llnl.gov/ 

5 http://stev.oapd.inaf.il/aesopus 
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Figure 1. Maps of Rosseland mean opacities showing the actual coverage of the tables in the log(T) - log(R) plane. While both plots refer to the same reference 
metallicity Z = 0.014, the left panel exemplifies the case of an "H-rich" table, the right panel corresponds to a "H-free" table, suitable to describe the opacity 
of a gas in which all helium has been burnt and converted into carbon and oxygen in the same proportions. 



while the helium content assumes Ny values (typically Ny = 10), 
ranging from Y = to Y = 1 - Z. In addition, for each Y value, 
we consider three combinations of C and O abundances, defined by 
the ratios R c = X C /(X C + X ) = 0.0,0.5, 1.0. All other elements 
are left unchanged. The "H-free" tables are specifically designed to 
describe the opacity in the He-burning regions. 

For both sets, each opacity table covers a rectangular region 
defined by the intervals 3.2 < log(r/K) < 8.7 and -8 < log(fl) < 1, 
where R = p/T^ in g cirr 3 , and = T/(10 6 K). To limit as much as 
possible the accuracy loss due to subsequent interpolation we adopt 
a fine grid spacing, with A log T = 0.01 for 3.2 < log(r/K) < 3.7, 
AlogT = 0.02 for 3.7 < log(r/K) < 8.7, and A logfi = 0.2. 
Two examples of opacity tables for both the "H-rich" and "H-free" 
cases are given in Fig. [T] The maps are obtained with the same 
grid spacing of the tables, which actually provides quite a dense 
and smooth description of the opacity all over the thermodynamic 
domain of interest. 

Interpolation over "H-rich" tables is perfomed in four dimen- 
sions, i.e. using R, T, X, and Z as the independent variables. While 
the interpolation in R and T is bilinear, we adopt a parabolic 
scheme for both X and Z interpolation. Interpolation over "H- 
free" tables is perfomed in five dimensions, i.e. involving R, T, 
Y, Rq = Xcl(Xc + Xo), and Z. Interpolation is bilinear in R and 
T, linear in R c , while we use as before a parabolic scheme for Z 
interpolation. 

Actually, bilinear interpolation in R and T does not preserve 
the continuity of the first derivatives, possibly introducing some in- 
stability into the convergence of the stellar model. For this reason, 
once the opacity tables are loaded, we compute and store the log- 
arithmic derivatives of the opacities (with respect to T and R) over 
the same grid of the opacity tables. Hence, the partial derivatives 
can be obtained by bilinear interpolation, in the same way as for 
the opacity values. While fulfilling the continuity requirements of 
the derivatives, this method has also the advantage of being quite 
fast. It should be also specified that, for each planned set of stellar 
tracks, the PARSEC code preliminarly loads a suitable number N z 
of opacity tables in order to follow in detail any significant change 
in the local metal content Z due, for instance, to the diffusion of 
heavy elements or dredge-up episodes. 



2.3 Equation of state 

For the equation of state (EOS) we make use of the FreeEOS code 
developed and updated over the years by A.W. Irwin, and freely 
available under the GPL licencq^. The FreeEOS package is fully 
implemented in our code and we may use it "on-the-fly", for differ- 
ent approximations and levels of accuracy. However, since the pre- 
tabulated version is sufficiently accurate for most of our purposes, 
we proceed by pre-computing suitable tables and by interpolating 
between them. 

The EOS calculation is performed accounting for the contribu- 
tions of several elements, namely: H, He, C, N, O, Ne, Na, Mg, Al, 
Si, P, S, CI, Ar, Ca, Ti, Cr, Mn, Fe, and Ni. For any specified distri- 
bution of heavy elements {Xj/Z}, we consider several values of the 
metallicity Z, and for each value of Z we pre-compute tables coin- 
taning all thermodynamic quantities of interest (e.g. mass density, 
mean molecular weight, entropy, specific heats and their deriva- 
tives, etc.) over suitably wide ranges of temperature and pressure. 
Exactly in the same fashion as for the opacity, we arrange a "H- 
rich" set containing N x = 10 tables each characterised by different 
H abundances, and a "H-free" set consisting of 31 tables, which 
are designed to describe He-burning regions. In practice we con- 
sider 10 values of the helium abundance, from Y = Oto Y = 1-Z. 
For each Y we compute three tables with C and O abundances de- 
termined by the ratios R c = X C /(X C + X ) = 0.0, 0.5, 1.0. 

Multi-dimension interpolations (in the variables Z, X or Y and 
Rc) are carried out with the same scheme adopted for the opaci- 
ties. All interesting derivatives are pre-computed and included in 
the EOS tables. 

Our procedure is to minimize the effects of interpolation by 
computing a set of EOS tables exactly with the initial metallicity 
and partition of the new set of tracks. This set is then inserted into 
the EOS database for interpolation when the global metallicity Z 
changes during the evolution. 

Before proceeding with the calculation of evolutionary tracks, 
we extensively tested the accuracy of the interpolation method. 
Fig. [2] illustrates the typical results from these tests: it compares 
evolutionary tracks computed with either the interpolation within 
the tabulated EOS, or a direct call to FreeEOS for each mesh point. 



6 http://freeeos.sourceforge.net/ 
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As can be appreciated, the differences in the HR diagram are clearly 
negligible. The same applies to other quantities in the tracks, such 
as the lifetimes and the luminosity at the RGB tip. 

2.4 Nuclear reaction rates 

Our nuclear network consists of the p-p chains, the CNO tri-cycle, 
the Ne-Na and Mg-Al chains, and the most important a-capture 
reactions, including the a-n reactions. The network solves for the 
abundances of N a = 26 chemical species: 'H, D, 3 He, 4 He, 7 Li, 
8 Be, 4 He, I2 C, I3 C, 14 N, 15 N, 1S N, 17 N, ls O, 19 F, 20 Ne, 21 Ne, 22 Ne, 
23 Na, 24 Mg, 25 Mg, 26 Mg, 26 AT, 26 AF, 27 Al, 28 Si. The latter nucleus 
acts as the "exit element", which terminates the network. In total 
we consider 42 reaction rates, listed in Tab. 12.41 These are the rec - 
ommended rates in the JINA reaclib database I Cvburt et alj|2010h . 
from which we also take the Q- value of each reaction. 

The electron s cree ning factors for all reac tions are those from 
iDewitt et al] i 19731) and lGraboske et all Jl973l) . The abundances of 
the various elements are evaluated wit h the aid of a s emi-im plicit 
extrapolation scheme, as described in iMarigo et al.l feOOlh . This 
technique, which does not assume nuclear equilibrium, is a con- 
venient compromise between the higher accuracy typical of the ex- 
plicit scheme, and the better stability of the solution guaranteed by 
the implicit scheme. 



2.5 Neutrino losses 



Energ y los ses by electron neutrinos are taken from lMunakata et al] 
dl985h and lltoh & Kohvam3 d 1 983t> . but for p lasma neutrinos, for 
which we use the fitting formulae provided by |Haftetal.l l ll994h . 



2.6 Convection 

2.6.1 Overshoot from the corrective core 

The extension of convective boundary of the core is estimated by 
means of an algorithm which t akes into account ov ershooting from 
the central convective region l lBressanetaT]|l98ll) . In our formal- 
ism the main parameter describing the overshoot is the mean free 
path of convective bubbles across the border of the convective re- 
gion, expressed in units of pressure scale heigh t, i.e. A c . Impor- 
tantly, this parameter in the lBressan et al] dl98ll) formalism is not 



equivalent to others found in literature. For instance, the overshoot- 
ing scale defined by A c = 0.5 in the Padova formalism roughly 
corresponds to the 0.25 pressure scale h eight above the con vec- 
tive border, adopted by the Geneva group dMevnet etal .11 19941 and 
references therein) to describe the same physical mechanism, i.e. 
A° = 0.25. 

It is well known that with this scheme it is not possible to use a 
unique choice of the overshoot parameter independently of the stel- 
lar mass. This is particularly true for main sequence (MS) stars in 
the transition between models with radiative and convective cores, 
i.e. in the mass range between M ~ 1.0 M Q and M ~ 1.5M Q , de- 
pending on the adopted chemical composition. In this mass range, 
a strict application of the above overshoot criterion give rise to the 
development of a convective core that grows too much with the evo- 
lution, producing a turn-off morpholog y which is not obser ved in 
the HR diagram (see the discussion in lAparicio et al ] |l99d) . This 
is true also for the solar model for which the convective core would 
persist up to the present age of 4.6 Gyr, which is not supported by 
current interpretation of helioseismology observations. 

We thus adopt a variable overshoot parameter A c to describe 
the development of overshoot, in the transition from main sequence 
models with radiative and convective cores: 

• A c is zero for stellar masses M < M 0[ , 

• In the range Moi < M < M02, we adopt a gradual increase of 
the overshoot efficiency with mass, from zero to A max . This because 
the calibration of the overshooting efficiency in this mass range 
is still very uncertain due to the scarcity of stellar data provided 
by the oldest interme diate-age and old open clusters. Some works 
dAparicio etai]|l990l) , however, indicate that this efficiency should 
be lower than in intermediate-mass stars. 

• For M > M 02 , we adopt A L . = A max . 



At variance with previous rele ases l lBressan et all 1 19931 : 

iGirardi et al]|2000l : iBertelli et al]|2008l) . where the limiting masses 
M i and M 02 were fixed (M i = 1 M G and M 02 = 1.5 M G ), we 
now adopt a transition region that depends on the initial chemical 
composition of the model set. This allows a better description of the 
effects of the elemental abundance on this phenomenon, which is 
needed because the onset of convective core overshoot leaves clear 
signatures that can be tested in observational Colour-Magnitude di- 
agrams (CMD). The transition masses are automatically defined by 
running a preliminary subset of stellar tracks without diffusion, to 
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Table 2. Nuclear reaction rates adopted in this work. 



Reaction 



Reference 



p(p,yS + v)D 

p (D , y) 3 He 

3 He( 3 He,y)2p+ 4 He 

4 He ( 3 He , y) 7 Be 

7 Be(e~ ,y)\A 

7 Li (p.yl'He + 4 He 

7 Be(p,y) 8 B 

12 C(p,r) 13 N 

13 C(p,y) 14 N 

14 N(p,y) 15 

15 N(p,y) 4 He+ 12 C 

15 N (p , y) 16 

16 0(p,y) 17 F 

17 0(p,y) 4 He+ 14 N 

17 0(p,y) 18 F 

18 0(p,y) 4 He+ 15 N 

18 0(p,y) 19 F 

19 F(p,y) 4 He+ 16 

19 F(p,y) 20 Ne 

20 Ne(p,y) 21 Na 

21 Ne(p,y) 22 Na 

22 Ne (p , y) 23 Na 

23 Na(p,y) 4 He + 20 Ne 

23 Na (p , y) 24 Mg 

24 Mg (p , y) 25 A1 

25 Mg (p , y) 26 A1S 

25 Mg (p , y) 26 Al m 

2£ Mg (p , y) 27 A1 

2£ AlS(p,y) 27 Si 

27 A1 (p , y) 4 He + 24 Mg 

27 A1 (p , y) 28 Si 

4 He(2 4 He,y) 12 C 

12 C( 4 He,y) 16 

14 N( 4 He,y) 18 F 

15 N ( 4 He , y) l9 F 

l6 ( 4 He , y) 20 Ne 

18 0( 4 He,y) 22 Ne 

20 Ne( 4 He,y) 24 Mg 

22 Ne( 4 He,y) 26 Mg 

24 Mg( 4 He,y) 28 Si 

l3 C( 4 He,n) 16 

17 O( 4 He,n) 20 Ne 

18 0( 4 He,n) 21 Ne 

21 Ne ( 4 He , n) 24 Mg 

22 Ne( 4 He,n) 25 Mg 

25 Mg( 4 He,n) 28 Si 



Cyburt et al. (20101 
Descouvemont et al. (2004) 
Angulo et al. (1999) 
Descouvemont et al. (2004) 
Caughlan & Fowler (1988) 
Descouvemont et al. (2004) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Imbriani et al. (2005) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Chafa et al. (2007) 
Chafa et al. (2007) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Iliadis et al. (2001) 
Hale et al. (2002) 
Hale et al. (2004) 
Hale et al. (2004) 
Iliadis etal. (2001) 
Iliadis etal. (2001) 
Iliadis etal. (2001) 
Iliadis et al. (2001) 
Iliadis et al. (2001) 
Iliadis etal. (2001) 
Iliadis etal. (2001) 
Fvnbo et al. (2005) 
Buchmann (1996) 
Gorres et al. (2000) 



Wilmes et al. (2002) 
Angulo et al. (1999) 
Dababneh et al. (2003) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Caughlan & Fowler (1988) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 
Angulo et al. (1999) 



check for the presence of a persistent convective unstable nucleus, 
without considering the overshoot (i.e., using the Schwarzschild 
criterium for convective borders). Since during the PMS phase the 
core may be already convective, we define the lower limit Moi to 
be the minimum mass that is able to maintain this unstable core 
after a 20 % of the initial hydrogen fraction has been burned in the 
center. 

Moi varies with both metallicity and helium content. The bot- 
tom panel of Fig. [3] illustrates the behaviour of M i as a function 
of metallicity for the main set of models to be described in this 
paper (see Sect. [5}. Notice that the resolution in the Moi determi- 
nation is of 0.05 M Q , which is the typical mass spacing we adopt 
for low-mass stars. The minimum values of Moi , of ~ 0.95 Msun, 
occur at the extremes of very low and very high metallicities. Moi 
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Figure 3. The behaviour of a few critical masses as a function of metal- 
licity, for our scaled-solar models following the Y = 0.2485 + 1.78 Z en- 
richment law (see Sect. |ffl. From bottom to top: the minimum mass that 
maintains a persistent convective core during H burning, Moi (see text for 
more details); the mass above which core overshoot is taken at the maxi- 
mum efficiency Mqi; the minimum mass for a model to ignite central He 
non degenerately, Mh c f for both the overshoot and the no-overshoot cases; 
finally, in the upper panel, the minimum mass of the stars that ignites C in 
a non electron degenerate core, Mim ■ The typical resolution in determining 
these mass limits is of 0.05 M Q for Moi and Mh c f, and of 0.2 M Q for Mim- 
The curves are polynomial fits to the corresponding values. 



is found to be 1.1 M Q at solar metallicities, and reaches a maximum 
of 1.15 M Q at Z = 0.02. These values are in good agreement with 
the lowest mass of a non-d iffusive model with a convective core, 
equal to 1.14 M Q , found bv lMichaud et al .1 d2004l) for assumed so- 
lar metallicities of eitherZ = 0.0175 orZ = 0.0199. The decreasing 
trend of Moi at lower metallicities can be explained by the larger 
luminosity at given mass, that causes the star to have a more con- 
centrated nucle ar energy source he nce a higher radiative gradient 
in the core (e.g. lMarigo et al.l200ll) . 

For the maximum overshooting efficiency we adopt A max = 
0.5, i.e. a moderate amount of overs hooting, which coin cides 
with the values ado pted in the previous iBertelli et al ] lfl99l) and 
iGirardi et all fcOOOh models. This corresponds to about 0.25 H P of 
overshoot region above the convective border found in other com- 
mon formalisms. 

M02 is always set to be Moi +0.3 M G . This choice is motivated 
by the modelling of the open cluster M 67 (see Sect.|7]and Fig.|16ll, 
which indicates overshooting as efficient as A c ^ 0.5 already at 
masses of ~ 1.3 M Q , for solar-metallicity stars. This choice is also 
supported by the SMC clu ster NGC 419 (Sect.|7]and IGirardi et all 
l2009l : lKamath et alfcOlfjh . in which the turn-off probes masses be- 
tween ~ 1.65 and 1.9 M Q . 

Furthermore, we assume that the overshooting efficiency A c 
increases linearly with mass between Moi and M02, in order to 
ensure a smooth enough transition between the properties of stars 
with radiative and convective cores. This prescription has to be con- 
sidered as a conservative approach, that will be revisited in future 
works (e.g. Rubele et al., in prep., and Rosenfield et al., in prep.). 
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The theoretical difficulties in defining the efficienc y of overshoot- 
ing in this transition region are well-known (e.g. lAparicio et al .1 
Il990l) . From the observational side, the recent indications from 
asteroseis mology are still ambiguo us: while the observations of 
a Cen A dde Meulenaer et al" suggest negligible overshoot- 

ing in solar-metallicity stars of mass ~ 1.1 M Q , recent astero- 
seismology studies of the nearby old low-mass star HD 203608 
dPeheuvels et al.ll201Clh . with [Z/X] a -0.5, indicate the action of 
overshooting (with a m = 0.17, which corresponds to a m ^ 0.32 
in our formalism) at masses as low as 0.95 M Q , which is probably 
just slightly above the Moi limit. Clearly, the behaviour of over- 
shooting in the transition region from M 0[ to M 02 deserves more 
detailed investigations. 

Another important difference with respect to previous releases 
of Padova evolutionary tracks is that the region of overshoot is con- 
sidered to be radiative. In our original scheme this region was as- 
sumed to be adiabatic while the majority of other authors took the 
overshoot region to be radiative. With the former assumption the 
model adjusts the internal radiative temperature gradient in such 
a way that the unstable region is smaller than in the case of a ra- 
diative overshoot region. Thus the effects of this new assumption 
are an overall larger mixed core. From a numerical point of view 
this assumption provides a faster convergence, because the original 
formalism causes some instability during the convergence process, 
that must be controlled. 

In the stages of core Helium burning (CHeB), the value A c = 
0.5 is used for all stellar masses. This amount of overshooting dra- 
matically reduces the extent of the br eathing pulses of c onvection 
found in the late phases of CHeB (see lChiosi et al.|[l992l) . 
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Figure 4. The luminosity both at the maximum (solid lines) and minimum 
("dip", dashed lines) of the bump feature along the first-ascent RGB, for 
a few low-mass tracks of (Z = 0.0005, Y = 0.249) computed with A c = 
0.05 H P (red) and A e = 0.5 H P (blue). In models with high A e the RGB 
bump is typically ~0.3 mag fainter than in model with negligible envelope 
overshooting. 



2.6.2 Overshoot in the convective envelope 

lAlongi etal] ( 1 199 ll) were the first to consider the possibility that 
the base of the convective envelope can give rise to a sizable over- 
shoot region. In the past, two important observational effects have 
been studied in relation with this phenomenon, namely the loca- 
tion of the RGB Bump in the red giant branch of low mass stars 
(globular clusters and old open clusters) and the extension of the 
blue loops of intermediate mass stars. Both effects have been found 
to be better explained by a moderate amount of overshoot, with a 
typical extension below the border of about 0.25 - l.OHp. Against 
this possibility it has always been argued that the calibration of 
the solar model do not require a sizable overshoot region because 
the transition between the fully adiabatic envelope and the radia- 
tive underlying region in our Sun, is already well reproduced by 
models without overshoot. This however does not exclude the pos- 
sibility that just below the fully adiabatic region convection may 
penetrate in form of radiative fingers that are even able to induce 
a significant mixing. Very recently it has been argued that a mech- 
anism of this kind could even provide a better agreement with the 
physical state of ma tter in this transition region derived from so- 
lar oscillations data l IChristensen-Dalsgaard et aT]|201 ll) . The size 
of this effect has been recently estimated to be A e ~ 0.4 H P but it 
is also consistent with a larger value A c ~ 0.6 H P which is in very 
good agreement with the one adopted since lAlongi et al.Ul99lh . As 
shown in Sect.[3]below , our own solar calibration suggests a value 
of A e (Adiabati c) < 0.05 H P . We have not yet inc luded the crite- 
rion adopted bv lChristensen-Dalsgaard et al .1 d201ll) . In view of the 
above results we will mai ntain in the present paper the the same 
prescr iptions adopted by lAlongi et al.1 jl99lh and iBressan et al.l 
jl993l) : A e = 0.05 in the envelope of stars with M < Moi and 
A c = 0.7 in the envelope of stars with M > Mo2- In between we 



adopt a smooth transition similar to that used for the central over- 
shoot in between these two limiting masses. Moreover, the temper- 
ature gradient in the envelope overshoot region is kept equal to the 
radiative one. 

It is worth recalling that, contrary to the overshoot in the cen- 
tral regions, the envelope overshoot has almost no effects on the 
evolutionary properties of the stars, though it may affect some ob- 
servable properties, as already discussed. In this respect, an im- 
portant observational evidence that we will investigate in the next 
future is the effect of the envelope overshoot on the abundance of 
light elements (e.g. 7 Li), during and after the PMS evolution. To 
this purpose one should also take into account that early main se- 
quen ce stars could still suff er from a non negligible mass accre- 
tion dDe Marchi et al.|[201 ll) and that envelope overshoot, during 
the P MS phase, could be m ore efficient than in the main sequence 
phase jMolaro et alj|2012h . However this is beyond the scope of 
this paper and will be addressed in a forthcoming investigation. 

Another important effect of envelope overshooting is that of 
lowering the luminosity level of the bump along the first-ascent 
evolution of low-mass stars. Fig. [4] compares the bump lumi- 
nosity between models computed with A e = 0.05 H P (red) and 
A c = 0.5 H P (blue) for a few low-mass tracks of composition 
(Z = 0.0005, Y = 0.249). In models with the higher A e , the RGB 
bump is typically ~ 0.3 mag fainter than in models with negligible 
envelope overshooting. By comparison, the observed RGB bump 
in globular clusters is about 0.2 to 0.4 mag fainte r than that pre- 
dicted by mo dels without envelo pe overshooting dPi Cecco et al.l 
l2010l : see also lCassisi et al.|[201ll) . though at the higher metallici- 
ties this result depends on the adopted metallicity scale. This sim- 
ple comparison favours larger values of A e in low-mass metal-poor 
stars, than here adopted based on the solar calibration. Needless to 
say, this argument will be the subject of further calibrations. 
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2.6.3 Breathing Pulses and Semiconvection 



2.8 Diffusion 



During the end of the core-He burning phase, models undergo a 
well known instability, the breathing pulses. In a series of pulses of 
the helium core a non negligible fraction of the helium in the sur- 
rounding stable regions is brought and burned into the core. This 
causes an increase of the helium burning lifetime which is not very 
significant by itself. But this has a dramatic effect on the duration 
of the subsequent AGB phase, which can be shortened by a sig- 
nificant factor. Given that this bears on the luminosity function of 
the brightest M and C stars, it is important to deal with this mech- 
anism in the more accurate and physically sound way as possible. 
At present there is no way to resolve the issue of whether this ef- 
fect is numerical or real, because it happens independently from 
the adopted mixing scheme, i.e. with mild overshoot or with semi- 
convection. Work is in progress to calibrate the efficiency of this 
mechanism with observations. 



2.6.4 Temperature gradient in the convective region 

The energy transport in the convecti ve regions is d e scribe d ac- 
cording to the mixing-length theory o f lBohm-Vitensel dl958h - The 
super-adiabatic gradient is maintained until its difference with re- 
spect to the adiabatic one decreases below VEicmcnt - V A di < 10~ 6 . 
While in the convective core this condition is always fulfilled, in 
the convective envelope it is fulfilled only near the bottom of the 
region. The mixing length parameter q-mlt is fixed by means of the 
solar model calibration performed below (see Section[3]l, and turns 
out to be q-mlt = 1-74. As already stated the temperature gradient 
in any overshoot region is assumed to be the radiative one. 



2.7 Mass meshing and other 

The stellar structure equations are solved maintaining the algorithm 
of the fitting method. The fitting mass is typically set at M/M Tol = 
1 - 10~ 5 , but it can be regulated on the basis of the corresponding 
temperature or depending on other needs. An adaptive mesh is used 
in the inner structure with a number of points that is about 1500 
during hydrogen burning, 3000 during He burning and 5000 during 
the thermally-pulsing AGB (TP-AGB). 

As for the atmosphere we adopt a plane-parallel grey model 
where the temperature stratification is given by a modified Edding- 
ton approximation for radiative transport: 



-Tt B [r + q(r)) 



(1) 



where r(r) is the optical depth, and q(r) is the Hopf function. As- 
suming the gas pressure P as the independent variable, the differ- 
ential equation dr/dlogP is solved by using a standard predictor- 
corrector method, starting from t = 0, where the gas pressure van- 
ishes and the total pressure is concides with the radiation pressure, 
to r = f , the latter being the optica l depth of the photospher e. The 
overall method is fully described in lKippenhahn et aT] l ll967l) . 

We explicitly follow the evolution of 26 elements, listed in 
Sect. 12.41 The network of nuc lear reactions is solved w ith a semi- 
implicit extrapolation scheme feader & Deuflhardl 19831) that treats 
all nuclear reactions together and withou t any assumption on nu- 
clear equilibria (see iMarigo et alj|200lh . At each time step the 
abundance equations are solved after each Henyey convergence, 
and particular care is paid both to mantain sufficiently small abun- 
dance variations (by adopting a suitably short time step), and to 
preserve the number of nucleons. 



Microscopic diff usion is included following the implementation by 
Salasnic iJdl999T) . The diffusion coefficients are calculated follow- 
ing lThouletal.1 dl994l) and he corresponding system of second 
order differential equations is solved together with the chemistry 
equation network, at the end of each equilibrium model. Diffusion 
is applied to all the elements considered in the code in the approxi- 
mations that they are all fully ionized. 

As already extensively discussed by IChabover etafl feOOll) . 
our knowledge of the diffusion process is still incomplete because 
if applied with the same diffusion constants used for the Sun in 
metal poor stars of similar mass, they lead to surface abundance 
that are in disagreement with observations. Indeed, on the basis 
of accurate abundance determi n ations in stars of glob ular cluster 
NGC 6397 l lGratton et al.ll200ll) . IChabover et aljj200lh concluded 
that microscopic diffusion should be fully inhibited in the external 
layers of metal poor stars at least down to a depth of about 0.005 M Q 
from the photosphere and partially inhibited down to 0.01 M a from 
the photosphere. 

In the case of metal rich stars, microscopic diffusion in the ex- 
ternal layers of stars of this masses is already inhibited by the more 
extended external convection. In general it is believed that, when 
one considers more massive stars with less extended surface con- 
vective regions, the effects of diffusion are negligible because of the 
much shorter stellar evolutionary times. Even in this case, however, 
it can be shown that if the mesh spacing in the external layers is 
kept suitably small for accuracy purposes, diffusion can noticeably 
change the surface composition even for masses well above those 
of globular clusters if the external convection disappears. This in 
turn changes the interior opacities and produces somewhat unre- 
alistic paths in the HR diagram, indicating that diffusion must be 
inhibited at the surface of these stars. 

Other criteria can be used to deal with the particular prob- 
lem of inhibiting molecular diffusion in those circumstances where 
its application would lead to evident discrepancies with observa- 
tions. Surface abundances in stars without extended surface con- 
vective regions are very se nsitive to the effects o f atomic diffusion 
and radiative accelerations dTurcotte et al.ll 19981) . However in gen- 
eral it is found that some sort of extra-mixing beyond the base 
of the external convective layers must be added in order to re- 
produce the observed abundances (e.g. iRicher etaL I l2000l) . This 
extra-mixing, of unknown origin but whose effect is that of mod- 
erating/inhibiting other diffusive processes, is parameterized as a 
turbulent diffusion with a coefficient that is cal ibrated on the ob- 
serve d surface abundances of old stars (see e.g. lVandenBerg et al.l 
I2012L and references therein). Unfortunately, when the calibrating 
observable is the surface Li abundance, the result s of different in- 
vestigations, based on differe nt stars, do not agree dMelendez et al.l 
l20ld : lNordlander et alj2012l) . Furthermore, this calibration is chal- 
lenged even more by the discovery that early ma in sequence stars 
still s uffer from a non negligible mass accretion jDe Marchi et ajj 
1201 ll) . This tail of accretion, preceded by a very efficient turbulent 
mixing during the pre-main sequence phase that completely destroy 
Li, could reshape ou r view of the surface evolution of this element 
dMolaroeta"I]|2012l) , and likely require another different calibra- 
tion of the inhibiting mechanism. In view of all the abov e reasons 
and s ince we neglect for the moment radiative levitation dVauclaij 
Il983l) . we switch off the diffusion process when the mass size of the 
external convection falls below a threshold mass that we assume to 
be AM c „ n v =0.5 % of the total stellar mass. Moreover diffusion is 
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Figure 5. Comparison between the turn-off age as a function of the stellar mass for models with and without diffusion (dark lines); with the relative difference 
plotted in red (see scale at the right axis). The right panel shows the same as a function of the turn-off luminosity. Differences between diffusive and non- 
diffusive models become negligibly small (less than 0.5 % in age) for masses above ~0.9 M Q . 
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Figure 6. Comparison in the HR diagram between evolutionary tracks com- 
puted with different thresholds for turning-off diffusion, for a set of tracks 
of low masses and for a low metallicity of (Z = 0.0006, Y = 0.249). Diffu- 
sion is switched off as soon as the mass size of the surface convection falls 
below the indicated fractional values of the total mass. For reference, the 
case without diffusion is also shown. See text for more details. 

not considered in stars that develop a persistent convective core, i.e. 
when core overshoot is taken into account. 

Figure [5] illustrates the impact of diffusion on the hydrogen 
burning lifetimes in low-mass metal poor evolutionary tracks with 
chemical composition (Z = 0.0006, Y = 0.249). From this figure 
we see that, with this criterion, the effect of diffusion becomes very 
small for masses above ~ 0.9 M G , which corresponds to ages of 
about 6 Gyr. At the 12 Gyr ages typical of globular clusters (with 
turn-off masses slightly below ~ 0.8 M Q ), the differences between 



diffusive and non-diffusive models become significant: for instance, 
for the same turn-off luminosity, ages derived from non-diffusive 
models would be about 8 % older than those derived from diffusive 
models. These d ifferences are perfectly in line w ith those found 
by other authors dCastellanT & deri'Innocentl 19991 , and references 
therein). 

To explore the sensitiveness of the results to different values of 
the threshold for turning-off diffusion, we have computed also the 
cases with AM conv =0. 1 % and 0.05 % for the same chemical com- 
position (Z = 0.0006, Y = 0.249). All the above cases are displayed 
in Fig. [6] Decreasing the threshold mass delays the quenching of 
diffusion and produces a larger shift of the turnoff region toward 
cooler temperatures. This effect is more pronounced in the mass 
range between 0.85 M G and 0.75 M Q . At lower masses surface con- 
vection is always well developed while at higher masses convec- 
tion rapidly disappears. However, we notice that in the latter case 
this result is also due to the inhibition of diffusion throughout the 
whole star, an assumption that could be improved following recent 
investigations that underline the i mportance of turbulent m ixing in 
the more external regions (see e.g. lVandenBerg et al.ll2012T ). In the 
light of these considerations we plan a revision of the tracks in this 
mass and metallicity range in the near future. 



3 CALIBRATION OF THE SOLAR MODEL 

The comparison with the solar model is a necessary step to check 
the quality of the input physics and to calibrate some parameters 
that cannot be derived from the theory. 

For this purpose we have computed a large grid of tracks of 
one solar mass, from the PMS phase to an age of 4.8 Gyr, and vary- 
ing the initial composition of the Sun, Zi„j ( j a i and ^initial! the mixing 
length parameter q-mlt and the extent of the adiabatic overshoot 
at the base of the convective envelope A c . To calibrate these free 
parameters, the models were compared with a set of solar data ob- 
tained from the literature which are summarized in Table[3] The ref- 
er ence solar data use d here come from MDI observation s described 
in iBasu et all fcoOuf) . More recent data can be found in lBasu et al.l 
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Figure 7. Solar model obtained using the tabulated Fre e-EOS. Left pane l: interior relative differences (in the sense Sun minus model) in the squared sound 
speed 5cJ/cj . Solar values are obtained from MD I dat a (Basu et al. 2000). Symbols are as follows. Red cross es: our model; blue dia monds ;:|BasuetalJ (2009) 
model differences with respect to MDI data using Grevesse & Sauvall U99a) solar abundan ces; black squares: ISerenelli et all 120091) using Grevesse & Sauval 
1 19981) solar abundances; green tria ngles: ISerenellietalTl2009ll using lAsplund et alj 120091) solar abundances. Right panel: interior relative differences for the 
density profile. In this case only the lBasu et al l 120091) model is shown for comparison. Z,, Yj and ML refer to the initial metallicity and helium abundance and 
to the MLT parameter of our model, respectively. The quoted age includes the PMS lifetime which is about 40 Myr, determined by the stage on the ZAMS 
where the total gravitational luminosity falls, and remains, below a few per cent. Values obtained from this calibration are shown in Table[5] 



Table 3. Solar calibration. 



Solar data 


Value 


error 


source 


LoaO^ergs" 1 ) 


3.846 


0.005 


Guenther et al. a 992) 


R e (10 10 cm) 


6.9598 


0.001 


Guenfher et al. (1992) 


7eff o (K) 


5778 


8 


from L & Rq 


Z Q 


0.01524 


0.0015 


Caffauetal. (2011) 


Ye 


0.2485 


0.0035 


Basu & Antia (2004) 


(Z/X) Q 


0.0207 


0.0015 


from Z Q & Y e 


«ADl/Ro 


0.713 


0.001 


Basu & Antia (1997) 


Padi 


0.1921 


0.0001 


Basu et al. (2009) 


Cs,ADi/10 7 cm/s 


2.2356 


0.0001 


Basu et al. (2009) 


Model 


Tab-EOS 


Fly-EOS 




Age + (Gyr) 


4.593 


4.622 




^initial 


0.01774 


0.01774* 




^initial 


0.28 


0.28* 




"MLT 


1.74 


1.74* 




Ae 


0.05 


0.05* 




L(10 33 ergs-') 


3.848 


3.841 




R(10 10 cm) 


6.9584 


6.96112 




Teff (K) 


5779 


5775 




Zs 


0.01597 


0.01595 




Ys 


0.24787 


0.24762 




(Z/X) s 


0.02169 


0.02166 




Radi/Ro 


0.7125 


0.7129 




Padi 


0.1887 


0.1881 




Cs,ADi/10 7 cm/s 


2.2359 


2.2364 





Distribution of heavy elements from Caffau et al] 1201 ll) . + Age includes 
the pre main sequence phase. 'Values taken from the calibration obtained 
with tabulated EOS. 



d2009h from BiSON experiment. However since we wish to com- 
pare our model with previous ones in literature, in particular con- 
sidering the effects of changing the present day surface solar com- 
position, we have found more convenient to use lBasu et al .1 J2OO0I) 
data, with respect to which there are several comparison at different 
chemical composition made in literature. Besides verifying the per- 
formance of our solar model another goal of the comparison is to 
obtain the mixing length parameter, o-mlt, that will be used to com- 
pute all the stellar evolutionary sets. For this reason the calibration 
has been obtained exactly with the same set-up used for the calcula- 
tions of the other tracks, i.e. with tabulated EOS and opacities and, 
of course, using microscopic diffusion. However in order to check 
the effects of interpolating the EOS we have also computed a solar 
model with the "on-the-fiy" version of FreeEOS. In this model we 
have adopted the parameters of the best fit obtained with the tabu- 
lated EOS, and we changed only the solar age in order to match as 
much as possible the solar data. The parameters of our best model 
are listed in the lower part of Table [3] 

In the left panel of Fig. [7] we plot the relative variation of 
the squared sound speed Sc^/c^ as a function of the fractional ra- 
dius inside the Sun (red crosses and solid line). The sense of the 
comparison is (Sun-model)/Su n in both panels. The Solar values 
were obtained from M DI data dBasu et afl l200d) . In the same fig- 
ure we also show: the lBasu et al] i 2009 | ) model differences with 
respect to MDI data, obtained with the Grevesse & Sauvall 1 19981) 
solar abundances (bl ue diamonds); the Serenel li et al.l d2009h dif- 
ference profile using iGrevesse & Sauvall J 1998h solar abundances 
(black squares) and using lAsplund et al.l ( 120090 solar abundances 



(green triangles). In the right panel we show the relative variation 
of the densit y profile (red cross es and solid line) and we compare it 
only with the lBasu et al j ( I2OO9T) difference profile (blue diamonds). 
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Notice that the quoted age includes the PMS lifetime, which when 
determined by the stage on the ZAMS where the total gravitational 
luminosity goes to zero, amounts to about 40 Myr. Inspection of Ta- 
bleland of Fig.|7]shows that our solar model performs fairly well, 
taking also into account th at lowering the solar metallicity from the 
iGrevesse & Sauval dl998h values, had a severe impact in the per- 
formance of the solar model. Particularly interesting are the small 
values of the <5cJ/c;r and density profiles in the central radiative re- 
gions of the Sun, indicating that good agreement is possible als o 
with abundances lower than those of IGrevesse & Sauvall Jl998h . 
Toward the central region our model suffer of a a slightly lower 
sound speed and a slightly lower density than shown by the solar 
data. A similar though less pronounced problem is also present in 
the compar ison of models wi th Solar data extracted from the Bison 
experiment lBasu et al.l d2009h . 

For the track computed with the "on-the-fly" FreeEOS we did 
not repeat the calibration process thus it does not reproduce the 
basic solar data of Table|3]as well as in the previous case. However 
we see that we may also obtain a fairly good agreement with solar 
data, by slightly increasing the age. 

From the initial values of the metallicity and Helium abun- 
dance of the Sun, JWtiai, ^initial, and adopting f or the primordial He 
abundance Y v = 0.2485 jKomatsu et alj|20lH) , we obtain also the 
helium-to-metals enrichment ratio, AK/AZ = 1.78. This value will 
be used in the following to determine the default value of helium 
content for any given global metallicity. 



4 INITIAL CHEMICAL COMPOSITION OF THE 
EVOLUTIONARY TRACKS 

As already specified in Section |2~T1 o ur selected reference distri - 
bution of heavy elements is taken from lGrevesse & Sauva 3i ll998h , 
except for a subset of species for which we adopt the recomm ended 
values according to the latest revision bv lCaffauetailfeOllh . 

For the present release, other distributions have been consid- 
ered such as those typical of the Large Magellanic Cloud, of mas- 
sive early-type galaxies and of a-depleted elements. 

The helium and metal mass fractions are chosen according to 
the relation Y = 0.2485 + 1 .78 Z obtained from the solar calibration. 
Sets with varying helium at fixed metallicity are also considered. 

We remind here that when we change the fractional abundance 
by number of heavy elements to obtain a new distribution, their 
fractional abundance by mass is re-normalized in such a way that 
the global metallicity, Z, is kept constant. 

Full sets of opacity and EOS tables are recomputed for each 
different chemical composition. Chemical compositions for which 
we have already computed full sets of evolutionary tracks are 
shown in Table [4] More sets are being computed for cr-enhanced 
distributions. 



5 STELLAR TRACKS 

A complete set of evolutionary tracks in the HR diagram is shown 
in Fig. [8] for the chemical composition Z = 0.008, Y = 0.263 with 
[a/Fe] = 0. Such plots are available for all sets listed in Table|4] 

5.1 Mass range 

In this paper we present the evolution of stars with initial mass be- 
tween 0.1 M G and 12 M e . Higher stellar masses will be considered 



Table 4. Initial chemical composition of the sets ready at the time of this 
writing 
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in a future paper. The mass spacing is smaller than in previous re- 
leases, being typically of 0.05 M Q in the range of low mass stars and 
~ 0.2 M G for intermediate-mass stars. This allows a good descrip- 
tion of the transition between low- and intermediate-mass stars, and 
between intermediate- and high-mass stars. 

5.2 Pre Main Sequence 

The PMS evolution is a new characteristic of this release and it is 
implemented as follows. A stellar model with suitable mass and ini- 
tial composition is obtained along the Hayashi track and artificially 
expanded until its central temperature falls below T c « 10 5 K. At 
this stage no nuclear reactions are active and the model is fully 
convective and homogeneously mixed. From this configuration, we 
allow the star to evolve at constant mass along the PMS contrac- 
tion and subsequent phases. The PMS phase of Z = 0.008, Y = 
0.263, [ff/Fe] = tracks is shown in Fig. [8] 

The typical PMS lifetimes are illustrated in Fig. [9] as a func- 
tion of stellar mass and for two different metallicities. It can be 
noticed that the behaviour as a function of mass is very regular 
just for masses higher than about 1 M Q . For smaller masses, the 
higher relative importance of the deuterium burning, and the onset 
of completely convective stars at M < 0.4 M G , cause the behaviour 
to be less monotonic. Also, it should be mentioned that the defi- 
nition of what is PMS phase is not so straightforward, so that in 
some cases we find some small glitches in the mass-PMS-lifetime 
relation, even at higher masses. These features seem to be less of a 
problem while making the isochrones for the PMS phase. 

5.3 Other evolutionary stages 

During the PMS evolution a number of low temperature reactions 
change the abundance of some light elements and the model is no 
longer homogeneous. The end of this phase is also characterized 
by a rapid reassessment of the internal structure during which the 
evolutionary track in the HR diagram shows characteristic patterns 
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Figure 8. Evolutionary tracks in the HR diagram for the selected chemical composition Z = 0.008, Y = 0.263 and with [ff/Fe] = 0. The upper panels 
show the tracks for intermediate- and high-mass stars (IMS+HMS) split into the PMS and later evolution; in the latter case the evolution is followed until the 
beginning of the TP-AGB, or up to C-ignition for M > 5.8 M Q . The low and very-low mass stars (LMS+VLMS) are shown in thebottom-left panels. In this 
case the evolution after the PMS is plotted either up to the tip of the RGB, or up to an age of 100 Gyr for M < 0.6 M Q . The bottom-right panel shows the 
complete suite of low-mass He-burning tracks computed for this set. The evolution starts at the ZAHB sequence and goes up to either the beginning of the 
TP-AGB for M > 0.6 M Q , or or up to the stages of post-early-AGB, and AGB-manque (and in some cases including the initial dimming towards the white 
dwarf sequence), typical of lower-mass HB stars that follow a hot-HB evolution (HHB). 



(loops) that depend on the initial mass and on the initial chemical 
composition. The above effects render difficult the identification of 
the Zero Age Main Sequence (ZAMS) point in an unambiguous 
way. To define the beginning of the MS we thus consider the point 
where the evolutionary speed in the HR diagram shows an abrupt 
drop, decreasing by more than two order of magnitudes in a very 



short time. This point is clearly detectable in all tracks and it will 
also be used to define the beginning of the MS when constructing 
the corresponding isochrones. 

Very low mass stars with mass M < 0.45 M Q are evolved 
up to an age of 100 Gyr. Some test models have been evolved up 
to the central Hydrogen exhaustion and subsequent core degener- 
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Mass (M Q ) 

Figure 9. PMS lifetimes, as a function of stellar mass and for two different 
metallicities. 

ate phase, when the model abandon the Red Giant Branch (RGB) 
because the Hydrogen shell reaches the surface of the star. These 
tracks show some interesting features but of course are not of any 
interest for single stellar evolution because the age where these fea- 
tures happens are by far larger than the Hubble time. 

Low mass stars, with mass between M ~ 0.5 M Q and an upper 
value (Mhsf) which mainly depend on the chemical composition 
and on the eventual extra mixing in the convective nucleus, ignite 
Helium in an electron degenerate hydrogen-exhausted core and un- 
dergo the so called He-flash. When the He-burning luminosity ex- 
ceeds the stellar surface luminosity by a factor of 200, the track 
is interrupted and the following evolution along the central Helium 
burning phase is computed separately. The He burning phase of low 
mass stars is re-started from a suitable Zero Age Horizontal Branch 
(ZAHB) model with the same core mass and internal and surface 
chemical composition as the last RGB model. The evolution is thus 
followed up to well developed thermal pulses on the AGB. The 
subsequent TP-AGB evolution up to the complete ejection of the 
envelope will be presented in an accompanying paper (Marigo et 
al., in preparation). 

In the construction of the initial ZAHB model we take into 
account that a fraction of the helium in the convective core has been 
burned into carbon during the He-flash. This fraction is estimated 
as the amount of nuclear fuel necessary to lift the core degeneracy, 
and is computed for each track. It closely follows the total mass of 
the core at the HeF, as illustrated in Fig. 1101 

For masses larger than M He F, we follow the evolution either 
until the first few (5-20) well developed thermal pulses along the 
AGB (in which case the star belongs to the range of intermediate- 
mass stars), or until central carbon ignition for higher mass models 
(in which case the star is broadly defined a massive star). For a 
given input physics (i.e. opacities, neutrino losses, etc.) the separa- 
tion mass between the two classes, M 1M , depends on the chemical 
composition and degree of convective core extra mixing. 

The values of Mh c f and Mm as a function of metallicity are 
illustrated in Fig. [3] for our grids of scaled-solar tracks. As can 
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Figure 10. Quantities that primarily determine the properties of He-burning 
low-mass stars, as a function of mass and for two solar-scaled chemical 
compositions. Bottom panel: The core mass at the moment of the He-flash. 
Middle panel: The same for the mass fraction of helium, AF C , that is burned 
into carbon to provide the energy to lift electron degeneracy in the core, 
between the RGB tip and ZAHB models. Top panel: The same for the 
mass fraction of helium in the envelope. 

be noticed, the behaviour of Mh c f closely follows the one of Moi 
(Sect. l27o~Tl l. with a maximum at intermediate metallicities. About 
the same happens for M m , but with the maximum values - of 
~ 6 M - moving to higher metallicities. 

5.4 a-enhanced mixtures 

Figure QT| exemplifies the effect of adopting tt-enhanced mixtures 
on the evolutionary tracks in the H-R diagram. There is a sys- 
tematic trend of the cc-enhanced tracks (blue line) to be somewhat 
warmer than the corresponding scaled-solar cases (red line), espe- 
cially along the RGB. 

In this respect, we call attention to the fact that the effect of the 
a-enhancement depends critically on the way the chemical mixture 
is effectively built. In general, given an [c/Fe] ratio, one has two 
options: a) either keeping the ratio [Fe/H] fixed while increasing 
the absolute abundances of the a elements, which leads to a net 
increase of the total metallicity Z; or b) keeping the metallicity Z 
fixed while depressing the abundances of the Fe-group elements 
and somewhat enhancing those of the a elements. The tracks pre- 
sented in Fig.[TT]are computed with the b) option. In the latter case 
the most relevant consequence is the effective increase of O, Mg, 
Ne, etc., while in the former case the most important effect is the 
depletion of the Fe-group elements. 

The changes in the abundances may produce important effects 
on evolutionary tracks, mainly due to opa city effects. We recom- 
mend to refer to iMarigo & Aringerl l l2009f) for a thorough discus- 
sion of this issue (in particular their sections 4.3, 4.4 and figures 
20, 21, 22). 

As a general rule we may expect that tt-enhanced tracks com- 
puted according to the a) option tend to be cooler than the corre- 
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Figure 11. Comparison between evolutionary tracks computed either with scaled-solar chemical composition , or cr-enhanced mixtures, while keeping the 
same initial total metallicity Z and helium content Y. Stellar masses (in M Q ) are indicated in each plot. 



sponding scaled-solar tracks (see e.g. IVandenBergetal.ll2012h be- 
cause of a net increase in the metallicity, while a-enhanced tracks 
computed according to the b) option tend to be cooler than the cor- 
responding scaled-solar tracks because the depression of the Fe- 
group elements implies a reduction of the H opacity, important 
absorption source at temperatures a 3000 - 6000 K. In fact the H 
opacity is highly sensitive to the number of free electrons, a signif- 
icant fraction of which is just provided by Fe-group elements. 



5.5 Comparison with previous releases 

In Fig.[l2]we compare our new tracks with the corresponding ones 
selected from a few previous releases from Padova, in order to give 
the reader an immediate impression of the similarities and differ- 
ences brought by the updates in input physics. The selected com- 
position is that typical of young populations in the Large Magel- 
lanic Cloud. The metallicity is Z = 0.008 for all the tracks but the 
Helium abundance Y may be slightly different among the different 
sets. 

Beside the presence of the PMS in the new tracks, a few gen- 
eral trends can be seen in the figure. 

Fist of all, the effect of diffusion is clearly visible in the tracks 
of lower mass. If we allow for a small shift in the ZAMS caused 
by the different input physics, we see that the new track runs cooler 
during the MS, and have a cooler turn-off. Another evident and 
important difference in low mass stars is that the base of the RGB 
is bluer in the new tracks. The two effects above will clearly affect 
the age determination of old populations, whenever the effective 
temperatures are used to const rain their ages . 

Moreover, with respect to lGirardi et alj d2000h . the RGB itself 
is less steep and runs at lower temperature s. In this phase the differ- 
ences with respect to lBertelli et al.1 d2008l) are less marked and the 
tracks have a very similar slope. 

In the intermediate mass stars (M = 5 M G ) the models show 
significant differences already from the MS. These differences are 
mainly due to the larger He abundance of the new tracks, and we 



limit the comparison only to the lBertelli et alj J2008h set, which has 
almost the same composition. In this case we can see that the new 
track runs slightly more luminous not only in the MS but also in 
the post-MS phases. This is mainly due to the different assumption 
about t he value of the real temperature gradient in the overshoot 
region. iBertelli et al.1 J2008h assume that the overshoot region is 
adiabatic. The consequence of this assumption in the equilibrium 
structure is that the Schwarzschild unstable core is less extended 
and thus it renders the extra mixing (with the same overshoot pa- 
rameter) less efficient. In other words assuming that the overshoot 
region is radiative is equivalent to assuming a slightly larger over- 
shoot parameter and, correspondingly, the star evolves at a higher 
luminosity because of the larger mixed core. 

We notice again that the red giant phase - in this case the early 
AGB - in the new tracks is significantly cooler than in lGirardi et al.1 
This difference is important because it bears on the mass- 
loss rates and thus on the subsequent evolution of the stars along 
the TP- AGB. 



5.6 Comparison with other databases 

There are several large databases of stellar evolutionary tracks 
available in the literature. All them differ in several important as- 
pects, including the adopted solar chemical composition and cali- 
bration, the sources of opacities and EOS, the coverage of evolu- 
tionary phases, the extension in the metal-helium content plane, etc. 
Comparing all these tracks is not among our goals. We limit our- 
selves to a few comparisons that point to the important similarities 
and differences between some recent libraries. 

Figure Q~3] compares the He core mass and the luminos- 
ity at the He-flash (or, equivalently, at the tip of the RGB), 
with those provided in the Dartmouth jDotter et alj 2008 ), Pisa 
dDell'Omodarme et alj|2012l) and BaSTI jPietrinferni et al]|2004h 
databases. It is evident that, although the initial helium content 
in not exactly the same among these tracks, the PARSEC mod- 
els agree very well with the most recent Pisa ones, and follow the 
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Figure 12. Comparison with previous releases of Padova ev olutionary tracks wit h Z = 0.008. Notice that He abun dances are sli ghtly differen t amon g the 
different sets illustrated: Y = 0.263 for PARSEC; Y = 0.25 for lFagotto et"ai] jl994h : Y = 0.25 for lGirardi et"aD <200fj : Y = 0.26 for lBertelli etal] l2008h . 




Figure 13. He core mass (left panel) and lumin osity (right panel) at the tip of t he red giant branch. The PA RSEC models for Z = 0.004, Y = 0.256 are 
compared with those fro m other authors, namely iDotter et aljf200l Danmouth), lDell'Omodarme et alj fimi Pisa), (both with Z = 0.004, Y = 0.250) and 
IPietrinferni eTai]i2004 BaSTI) (with Z = 0.004, Y = 0.251). 



same trend as the Dartmouth ones - with just small offsets, of a few 
0.001 Mq in the He core mass, and a few 0.01 dex in luminosity. 
This general agreement is very encouraging. It is also interesting 
to note that, although presenting offsets of similar magnitude, the 
BaSTI tracks have a different behaviour, with slightly smaller core 
masses being followed by slightly larger luminosities at the tip of 
the RGB. 

As an illustration of the possible differences in the 
HR diagrams. Fig. [T4l compares PARSEC with BaSTI tracks 
jPietrinferni et al.ll2004l) . using always the same chemical compo- 
sition of Z = 0.0198, Y = 0.2 734. This is the init i al sola r compo- 
sition in BaSTI, obtained from lGrevesse & Sauvall Jl998h and their 



calibration of the solar model, which results in a mixing length pa- 
rameter of oimlt = 1-913. The straight computation of PARSEC 
stellar evolutionary tracks with the same Z, Y, and q-mlt results 
in tracks very similar to BaSTI, as can be appreciated compar- 
ing the red (BaSTI) with the blue lines (PARSEC). However, the 
bulk of present PARSEC tracks are calculated with our own solar- 
calibrated «mlt = 1.74. PARSEC tracks with q-mlt = 1-74 are pre- 
sented in green. They appear significantly cooler than the BaSTI 
tracks, especially at the RGB. This significant difference results 
mainly from the lower metallicities assumed for the solar model 
in PARSEC, and the calibration of q-mlt- 

This example illustrates how risky is the straight comparison 
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Figure 14. Comparison with BaSTI tracks jPietrinferni et al.l :2004 ) in the HR diagram. The two panels show trac ks for masses of 1 .0 and 2 .6 M Q obtained from 
the BaSTI (red lines) database for the same fixed solar-scaled abundance of Z = 0.0198, Y = 0.2734 (based on iGrevesse & Sauvallll998h and q-mlt = 1.913. 
The blue lines are PARSEC scaled-solar tracks obtained for the same composition and ffMLT value as in BaSTI, whereas the green lines are the same but using 
the PARSEC calibration of the mixing length parameter, i.e. «mlt = 1.74. 



between tracks computed with a different solar composition and 
calibration. 



5.7 Data tables 

The data tables for the present evolutionary tracks are available in 
electronic format only through our websites at OAPE0. For each 
evolutionary track, the corresponding data file presents at least 50 
columns with the run of many quantities of interest along the track, 
as described in a suitable ReadMe file. 



6 ISOCHRONES 

For the sake of completeness we present here, together with the set 
of evolutionary tracks, also the corresponding isochrones. These 
isochrones do not contain neither the TP-AGB nor the post-AGB 
phases. These are discussed in accompanying papers dealing with 
the extension of the tracks into these advanced phases (e.g. Marigo 
et al., in preparation), and their calibration (basically in terms of 
third dredge-up and mass loss efficiencies) on the base of a large 
set of observational data. Although we can ensure a quick initial 
ca libration of the TP - AGB phase following the procedure defined 
by|M artgo &Girardilj2007l) . a thoroughful calibration may require 
a significant time before completion. We thus focus here our at- 
tention on the less advanced phases of stellar evolution which are, 
nevertheless, important thanks to the huge amount of observational 

data they can h elp to interpret . 

Following iBertelli et alj jl994h . isochrones are computed by 
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first dividing the evolutionary tracks in suitable homogeneous evo- 
lutionary phases. The isochrones are then constructed by interpo- 
lating points along missing stellar tracks, between the correspond- 
ing phases of the calculated grid masses. The points that separate 
the selected evolutionary phases along the tracks are automatically 
recognized by a suitable algorithm that takes into account both the 
morphology in the HR diagram and the interior physics. We store 
the full track information for the purpose of constructing isochrones 
even in non-conventional diagrams, like for example in the surface 
composition of some particular elements. 

Examples of isochrones in the HR diagram are shown in 
Fig. [15] for the sets with scaled-solar chemical composition that 
follow the relation Z = 0.2485 + 1.78 Z. For these sets, isochrones 
can already be built for any age from ~ 1 Myr to at least 20 Gyr 
(limited to the mass interval 0.1 < M < 12 M B ), and for any 
metallicity in the range 0.0005 < Z < 0.07 (corresponding to 
-1.49 < [M/H] < +0.78). We also note the presence of the PMS in 
the isochrones with the youngest ages. Mass loss by stellar winds 
has been considered only d uring the RGB of low-mass stars, us- 
ing the empirical formula bv lReimer si d 19751) multiplied by an effi- 
ciency factor j]. The latter is set as 77 = 0.2, w hich is representativ e 
of the modest values recently determined bv lMiglio etai] ( l2012h . 
based on the asteroseismic data for two star clusters in the Kepler 
fields. 

The initial point of each isochrone is the 0. 1 M Q model in the 
lower main sequence. The terminal stage of the isochrones is either 
the beginning of the TP-AGB for low- and intermediate-mass stars 
(M < Mim, ages > 10 8 yr), or C-ignition for more massive stars. 

Theoretical luminosities and effective temperatures along the 
isochrones are translated to magnitudes and colors using exten- 
sive tabulations of bolometric corrections and colors, as detailed 



in lGirardi et ail d2002l).lMarigo et alj l2008h . iGirardi et alj J2008h , 



and Rubele et al. (2012). The tabulations were obtained from con- 
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Figure 15. Examples of theoretical isochrones in the HR diagram. Left panel: A sequence of solar-metallicity isochrones for ages going from log(7/yr) = 6 
to 10.2 at equally spaced intervals of A log; = 0.1. Right panel: Three sequences of isochrones for a fixed age and at varying metallicty. The ages are 
log(f/yr) = 6.3, 8.0 and 10.1. The sequence of metallicities goes from [M/H] = -1.45 to +0.75 at equally spaced intervals of A[M/H] = 0.2. In all cases, 
the mass-loss parameter is assumed to be r\ = 0.2, and the sequences are completed down to 0.1 M Q . Notice the presence of the PMS phase in the youngest 
isochrones. 



volving the spectral energy distributions cont a ined in the libraries 
of ste llar s pectra of ICastelli & Kuruczl <2003l) . lAllard et alj feOOfj. 
l200ll) . and bluks et alj Jl994) with the response function of several 
medium- and broad-band filters. The refe rence spectrum for Vega- 
mag systems is taken from lBohlinl J2007l> . 

The bulk of stars modelled in optical filters, for r c ff be- 
twe en 50000 and ~ 400 K, have their colors based on 
the ICastelli & Kuruczl J2003h ATLAS9 models. Casagrande et al. 

have recently and accurately revised the temperature scale 
of nearby stars via the infrared flux method. It is extremely reassur- 
ing to note t hat, for the Sun, the color d ifferences between their cal- 
ibration and lCastelli & Kuruczl J2003h models are in general only 
=* 0.02 mag (which is also comparable to the random errors in the 
calibration; see their table 7), all the way from the B to the H band. 

The isochrones are initially made available through the CMD 
web interface at OAPE0. It allows isochrones to be retrieved for 
any given choice of age and initial metallicity, within the range of 
the calculated tracks. Output tables include the most relevant stellar 
parameters along the isochrones (e.g. L, T cff , g, radius, evolutionary 
phase), and the absolute magnitudes in the photometric system of 
choice. Isochrones are now provided in more than 30 photometric 
systems corresponding to major astronomical facilities and surveys. 



http://stev.oapd.inaf.it/cmd 



However, the list of photometric systems is being continuously ex- 
panded, and updated in response to changes into filter transmission 
curves and zero-points. We anticipate that the same web interface is 
being restructured in order to deal with the entire variety of stellar 
parameters stored in the PARSEC evolutionary tracks. 



7 DISCUSSION AND CONCLUDING REMARKS 

We have discussed the main features of the new code PARSEC that 
will be used to update and extend the Padova database of stellar 
evolutionary tracks and isochrones. In PARSEC we have modified 
and updated all the major input physics, including the equation of 
state, the opacities, the nuclear reaction rates and nuclear network, 
and the inclusion of microscopic diffusion. 

With the new code we have computed several sets of stellar 
evolutionary tracks that differ from the previous ones in a few ad- 
ditional aspec ts. (1) They a r e base d on a new reference solar com- 
position from ICaffau et al.1 d201 lh . The present code however al- 
lows the quick computation of tracks with modified metal distribu- 
tions (e.g. a enhancement or depletion), fully taking into account 
the changes implied to tables of opacities and equation of state. (2) 
We consider the evolution along the PMS phase. (3) A few changes 
have been introduced on the way the convection is treated. From 
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Figure 16. PARSEC isochrones overlaid in the CMD of the intermediate- 
age star cluster M 67. The distance modulus and colour excess were fixed 
as V - M v = 9.75 and Eb-v = 0.03, as shown in Fig. [JjJ next. The BV 
photometry is taken from Montgome ry et al] il993l black points). Red dots 
mark the "high-probability single stars" identified bv lSan dquist (20(3). The 
isochrones are for Z = 0.014, and for several values of overshooting param- 
eter A c . The ages are selected so as to reasonably fit the position of the 
subgiant branch, and especially its lower boundary which coincides with 
the locus of bona-fide single stars. Notice that the isochrones with A c < 0.4 
and A c > 0.6 fail in reproducing the position and shape of the turn-off re- 
gion. 

these tracks we obtain preliminary isochrones up to the beginning 
of the TP-AGB. 

Examples of the quality of the new isochrones are shown in 
Figs. [16] and [TTJ where we present preliminary fits of the colour 
magnitude diagrams of the Galactic open cluster M 67 and the SMC 
cluster NGC 419. 

M67 is the best example of an intermediate-age cluster in 
which the turn-off region presents the clear signatures of well- 
developed convective cores, in particular the sizeable gap in the 
concentration of stars at the termination of the main sequence, 
which is followed by the rapid contraction of the previously con- 
vective cores, soon after H-exhaustion at the centre. This cluster 
indicated the need for a quick increase of overshooting efficiency 
with mass, as advanced in Sect. 12.6.71 Indeed, Fig. 1 161 shows a se- 
ries of isochrones obtained from PARSEC tracks with overshooting 
parameter going from A c =0.1 to 0.7, for the near-solar composi- 
tion of (Z = 0.014, Y = 0.273) ([M/H] = -0.02). In all cases, the 
isochrone age was fixed so as to produce about the same descrip- 
tion of the subgiant phase, along the line that joins the bona-finde 
single stars located between the reddest termination of the turn-off 
region and the lower extremity of the fisrt-ascent RGB. The dis- 
tance and reddening of the isochrones is fixed so as to satisfactorily 
reproduce the mean position of the RGB, the lower main sequence, 
and the red clump of He-burning stars, as can be appreciated in the 
left panel of Fig. llifl . Therefore, all these isochrones differ mainly 

9 We recall that for the typical masses in these isochrones, the red clump 



in their description of the turn-off region. The figure reveals that 
isochrones with A c < 0.4 present a too blue turn-off, compared to 
the locus of bona-fide single stars in M67. Conversely, isochrones 
with A c > 0.5 have clearly too red a turn-off. This comparison sup- 
ports our initial choice for the value of overshooting efficiency for 
masses M > M 02 , as described in Sect. |2.6.T| More detailed work is 
underway, considering the full range of parameters (distances, red- 
dening, metallicity, binary frequency, etc.), as well as the detailed 
predictions about the occupation probability of different sections of 
the isochrones. 



We also note that M67 has been object of many recent 
analyses aimed at deriving accurate determinations of its age, 
and constraints on its chemical composition and convective effi- 



ciencv (see e.a.lMichaud et alj|2004; 


VandenBerg & Stetsonll2004|; 


VandenBera et al. 2007; Magic et al. 


l201fJ: iMowlavi et alj|2012h. 



These works were intentionally focused in the turn-off region of 
the CMD - and especially on its hook-like feature - and did not 
attempt to fit the RGB and red clump altogether with the same 
isochrone, as we do in Fig. [TTJ Our choice of a "global isochrone 
fitting" is based on a very simple consideration: our evolutionary 
tracks are mostly aimed at improving evolutionary population syn- 
thesis of resolved and unresolved galaxies, for which most observa- 
tions include very well sampled giant branches. Therefore, a good 
reproduction of the evolved phases in the CMD is as important in 
our case, as is the good reproduction of the turn-off region in those 
more focused works, where the main target is the calibration of 
physical processes and abundances inside stars. 

Another interesting object for a preliminary check of the 
efficiency of overshooting is the intermediate-age SMC cluster 
NGC 419. Its CMD (right panel in Fig. [J7J> presents clear evi- 
dences of multiple populations , namely a broad turn-off and a dual 
red clump l lGirardi et alj2009l) . The turn-off masses are around the 
M He F limit between intermediate- and low-mass stars, as evidenced 
by the double structure of the red clump. For a cluster of this kind, 
reproducing the shape of both the turn-off and red clump regions 
is possible only if the isoch rones have about the right amount of 
overshooting, as discussed in iGirardi et aljj2009l) . As illustrated in 
Fig-OS the default PARSEC isochrones (with A c > 0.5 in the mass 
range of interest for NGC 419) with ages spanning from 1.4 to 2.0 
Gyr are able to reproduce both the "golf-club" shape of the turn-off 
region, and the dual structure of the red clump - where the sec- 
ondary, fainter red clump is explained by the youngest isochrone. 

For the moment, our goal is just to call attention to the cor- 
rect reproduction of the main CMD features, including the shape 
of the turn-off(s), the subgiant branch, the RGB slope, and the red 
clump(s) in these two clusters. We anticipate that a detailed quanti- 
tative analysis of these and other similar clusters is being performed 
to calibrate the efficiency of overshoot in this mass interval (Rubele 
et al., in prep.). Younger stellar populations will be considered by 
Rosenfield et al. (in prep.). The evolution along the TP-AGB, and 
the corresponding isochrones, will be presented in an accompany- 
ing paper (Marigo et al., in prep.). 
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Figure 17. PARSEC isochrones overlaid in the CMDs of two intermediate-age star clusters. In both cases, the distance modulus and reddening were fixed "by 
eye" so as to reproduce the position of the red clump. Left panel: The Galactic open cluster M 67. The data are the same as in Fig.[l6] The isochrones are for 
Z = 0.014 and log(f/yr) = 9.6, distance modulus of V - M v = 9.75 and color excess = 0.03. Right panel: The core of the SMC cluster NGC 419. The 
photometry comes from archival HST images obtaine d with the Hi gh Resolution Channel of the Advanced Camera for Surveys onboard HST (GO-10396, 
PI: J.S. Gallagher), and reduced as described in iGirardi et all 120091) . In order to reduce the presence of blends and spurious objects, we plot only those stars 
with a sharpness of (sharpp 555w + sharp Fgl4w ) 05 < 0.1. The isochrones are for Z = 0.004, ages log(f/yr) = 9.15, 9.20, 9.25, 9.30 (from blue to red), distance 
modulus of F814W - M^i4w = 18.96 and color excess Efsssw-J'SIW = 0.09. 
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